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ABSTRACT 



The Alzheimer's Disease Neuroimaging Initiative (ADNI) recently added diffusion tensor imaging (DTI), among 
several other new imaging modalities, in an effort to identify sensitive biomarkers of Alzheimer's disease (AD). 
While anatomical MRJ is the main structural neuroimaging method used in most AD studies and clinical trials, 
DTI is sensitive to microscopic white matter ( WM) changes not detectable with standard MR], offering additional 
markers of neurodegeneration. Prior DTI studies of AD report lower fractional anisotropy (FA), and increased 
mean, axial, and radial diffusivity (MD, AxD, RD) throughout WM. Here we assessed which DTI measures may 
best identify differences among AD, mild cognitive impairment (MCI), and cognitively healthy elderly control 
(NC) groups, in region of interest (ROI) and voxel-based analyses of 155 ADNI participants (mean age: 73.5 ± 
7.4; 90 M/65 F; 44 NC, 88 MCI, 23 AD). Both VBA and ROI analyses revealed widespread group differences in 
FA and all diffusivity measures. DTI maps were strongly correlated with widely-used clinical ratings (MMSE, 
CDR-sob, and ADAS-cog). When effect sizes were ranked, FA analyses were least sensitive for picking up group 
differences. Diffusivity measures could detect more subtle MCI differences, where FA could not. ROIs showing 
strongest group differentiation (lowest p-values) included tracts that pass through the temporal lobe, and poste- 
rior brain regions. The left hippocampal component of the cingulum showed consistently high effect sizes for 
distinguishing groups, across all diffusivity and anisotropy measures, and in correlations with cognitive scores. 

© 2013 The Authors. Published by Elsevier Inc. All rights reserved. 



1. Introduction 

Alzheimer's disease (AD) is the most common type of dementia, af- 
fecting 1 in 8 people over age 65 in the U.S. alone (Alzheimer's Disease 
Association, 2012). Its prevalence is predicted to more than double in 
the next 40 years (Hebert et al., 2003). It is important to identify 
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individuals most likely to develop AD, so that those at greater risk can 
be treated earlier. One high-risk group consists of people with mild cog- 
nitive impairment (MCI) — a transitional stage between normal aging 
and AD. People with MCI convert to AD at a rate of about 10-15% per 
year (Petersen et al., 2001; Bruscoli and Lovestone, 2004). In addition 
to the more widely-accepted measures from anatomical MRI, PET, and 
CSF measures of pathology, one major neuroimaging study of AD - the 
Alzheimer's Disease Neuroimaging Initiative (ADNI) - recently incorpo- 
rated additional neuroimaging techniques including diffusion tensor 
imaging (DTI) (Jack et al., 2010; Jahanshad et al, 2010a; Zhan et al., in 
press). DTI is a variant of MRI that measures the diffusion of water mol- 
ecules in brain tissue. Here we set out to assess which standard DTI 
measures may best identify neuroanatomical differences between AD, 
MCI, and normal aging. In the end, DTI offers a range of measures that 
might be sensitive to pathology, including measures of brain connectiv- 
ity (Daianu et al., 2012, 201 3a,b; Nir et al., 2012; Prasad et al., 2013; Toga 
and Thompson, 2013). For these initial analyses, however, we aimed to 
analyze more traditional measures and maps that are perhaps most 
likely to be used in standardized multi-site DTI analyses, at least in the 
near future (Jahanshad et al., 2013). 



2213-1582/$ - see front matter © 2013 The Authors. Published by Elsevier Inc. All rights reserved. 
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MRl-based image analysis methods have long been used to track 
structural atrophy of the aging brain. MRI studies of AD reveal wide- 
spread neuronal loss and atrophy in the brain's gray matter, especially 
in medial temporal and hippocampal regions (Atiya et al., 2003; 
Chetelat and Baron, 2003; Thompson et al., 2003; Anderson et al., 
2005; Apostolova and Thompson, 2008; Bakkour et al., 2009; Risacher 
et al., 2009; Apostolova et al., 2010; Desikan et al., 2010a; Desikan 
et al., 2010b; Chiang et al., 2011; Weiner et al., 2012; Leung et al., 
2013). Beta-amyloid and tau proteins accumulate in the brain, leading 
to inflammation, neuronal atrophy and cell death (Braak and Braak, 
1991 ; Braak and Braak, 1995). As neurons are lost, white matter volume 
is also reduced, due to both myelin degeneration and axon loss in neural 
fiber tracts (Braak and Braak, 1996; Bartzokis, 2011; Braskie et al., 2011; 
Hua et al., 2013). Standard anatomical MRI is still the imaging technique 
most often used in AD studies and clinical trials, but DTI is sensitive to 
microscopic changes in white matter (WM) integrity not always 
detectable with standard anatomical MRI (Xie et al., 2006; Canu 
et al., 2010). Although this is debatable until more evidence is col- 
lected, some DTI changes may even precede and predict volume loss 
(Hugenschmidt et al., 2008; Nir et al., 2012), making it a potentially 
beneficial tool for capturing additional or complementary markers of 
early neurodegeneration. Carriers of some AD risk genes show differ- 
ences on DTI as young adults, decades before the typical age of onset 
ofAD(Braskie et al., 2011). 

Fractional anisotropy (FA) is perhaps the most widely accepted DTI 
measure and reflects how directionally constrained the diffusion of 
water is along axons. While higher FA values may imply more coherent 
or intact axons, or a higher degree of myelination, lower FA may reflect 
loss of WM integrity and injury. These physiological correlates of the DTI 
signal are widely accepted, but the differences may have other interpre- 
tations, especially where fibers cross (Leow et al., 2009; Zhan et al., 
2009). Mean diffusivity (MD) captures the average rate of diffusion in 
all directions, and generally increases with WM injury, especially if nor- 
mal barriers to diffusion are damaged (such as myelin sheaths on 
axons). Axial diffusivity (AxD) captures diffusion parallel to axonal fi- 
bers, while radial diffusivity (RD) reflects perpendicular diffusion. 
These measures are linked to axonal injury and demyelination, respec- 
tively (Song et al., 2003 ; Song et al., 2005). To date, numerous DTI studies 
of AD and MCI find that greater cognitive impairment, or poorer diagno- 
sis, is associated with lower FA in the corpus callosum, fornix, cingulum, 
superior longitudinal fasciculus, and inferior longitudinal fasciculus 
(Ukmar et al., 2008; Strieker et al., 2009; Mielke et al., 2009; Liu et al., 
201 1 ) and DTI measures correlate with widely used clinical or cognitive 
ratings including the mini-mental state exam (MMSE) (Bozzali et al., 
2002). 

Despite growing diffusion imaging evidence of AD-related WM 
changes, it is not clear which regions and DTI measures are the most sen- 
sitive for detecting diagnostic differences. In order to evaluate the power 
of drag trial treatment to counteract degeneration, optimizing statistical 
power for discerning differences and changes is crucial. We focused this 
current paper on cross-sectional differences in patients and controls, as 
there are a number of DTI measures, regions, and approaches that need 
to be compared and ranked in terms of their effect sizes for picking up 
group differences. We set out to rank the effect sizes for different DTI- 
based scalar measures in detecting differences in both white matter 
voxel-based analyses (VBA) and within regions of interest (ROIs). We 
first examined differences in DTI anisotropy and diffusivity measures be- 
tween groups of cognitively healthy normal elderly controls (NC), MCI, 
and AD patients in both voxel-based and ROI analyses. We also exam- 
ined the association of anisotropy and diffusivity maps with widely 
used clinical or cognitive ratings including the MMSE (Folstein et al., 
1975), the "sum-of-boxes" clinical dementia rating (CDR-sob) (Berg, 
1988), and the Alzheimer's Disease Assessment Scale-Cognitive (ADAS- 
cog) (Rosen et al., 1984). Finally, in a supplementary test, we compared 
our ROI results to ROIs extracted along the skeleton from the widely used 
tract-based spatial statistics (TBSS) method (Smith et al., 2006). Despite 



the popularity of FA, we hypothesized that we would find the highest ef- 
fect size and discriminative power for MD measures, as recently sug- 
gested in a review of DTI studies of AD by Clerx et al. (Clerx et al., 
2012). We also hypothesized that we would find the greatest differences 
in temporal lobe WM and the corpus callosum (CC), as the temporal lobe 
is usually the earliest region to be affected by amyloid and tau pathology 
in AD, and DTI studies are often better powered to find group differences 
in regions such as the CC where fiber coherence is highest. 

2. Materials and methods 

2.1. Clinical sample and demographics 

Baseline MRI, DTI, clinical, and neuropsychological data were 
downloaded from the ADNI database (http://adni.loni.ucla.edu). When 
the analysis was performed (September 2012), data collection for the 
ADNI2 project was still in progress. Here we performed an initial analy- 
sis of 155 participants from 14 data acquisition sites, of whom 44 were 
normal controls (NC), 88 amnestic MCI subjects, and 23 AD patients 
(see Inline Supplementary Table SI for distribution of subjects across 
sites). Unlike ADNI1, ADNI2 MCI participants include the enrollment 
of a new early MCI cohort (e-MCI; n = 62), with milder episodic mem- 
ory impairment than the MCI group of ADNI1. The MCI group of ADNI1 
is now referred to as late MCI (1-MC1; n = 26) in ADNI2. Levels of MCI 
(early or late) were determined using the Wechsler Memory Scale — 
Logical Memory II (Wechsler, 1987). We evaluated the 1-MCI and e- 
MCI groups both separately and as one large MCI group. Detailed inclu- 
sion and exclusion criteria are found in the ADNI2 protocol (http://adni- 
info.org/Scientists/Pdfs/ADNI2_Protocol_FINAL_20100917.pdf). 

Inline Supplementary Table SI can be found online at http://dx.doi. 
org/10.1 01 6/j.nicl.201 3.07.006. 

Each subject underwent cognitive evaluations. The Mini-Mental 
State Examination (MMSE) was used to provide a global measure of 
cognitive status, based on evaluating cognitive domains including ori- 
entation to place, orientation to time, registration, attention and con- 
centration, recall, language, and visual construction (Folstein et al., 
1975). The total score ranges from 1 to 30, with lower scores indicating 
impairment. The Clinical Dementia Rating (CDR) was also used as a 
global measure of dementia severity (Berg, 1988). The "sum-of-boxes" 
CDR (CDR-sob) score is the sum of 6 measures each assessing the de- 
gree of impairment in memory, orientation, judgment and problem 
solving, community affairs, home and hobbies, and personal care. The 
CDR-sob score ranges from 0 to 18 (no dementia to severe dementia, re- 
spectively). Finally, the Alzheimer's Disease Assessment Scale-Cognitive 
(ADAS-cog), a global measure encompassing memory, reasoning, lan- 
guage, orientation, ideational, praxis and constructional praxis (Rosen 
et al., 1984), was collected where scores range from 0 to 70 (no demen- 
tia to severe dementia respectively). In post-hoc analyses, we further 
homed in on specific cognitive domains using the available ADNI 
composite scores for executive function (ADNI-EF) (Gibbons et al., 
2012) and memory (ADNI-MEM) (Crane et al., 2012) derived using 
data from the ADNI neuropsychological battery. Detailed psychometric 
calculation protocols are available for download at https://ida.loni.ucla. 
edu/. ADNI-EF was calculated using a combination of WAIS-R Digit 
Symbol Substitution, Digit Span Backwards, Trails A and B, Category Flu- 
ency, and Clock Drawing scores (Gibbons et al., 2012), and ADNI-MEM 
was calculated as a composite of the Rey Auditory Verbal Learning 
Test (RAVLT), ADAS-Cog, and Logical Memory data (Crane et al., 2012). 

Demographics and diagnostic information for the participants are 
shown in Table 1 . Diagnostic groups did not differ in age, however, edu- 
cation, an AD risk factor (Sattler et al., 2012), was marginally significant 
between controls and AD. As would be expected, clinical measures 
that index cognitive decline (MMSE, ADAS-cog, CDR-sob, ADNI-MEM, 
ADNI-EF) did show significant graded differences between groups. 

We further assessed whether these measures revealed more fine- 
grained differences between the 1-MCI and e-MCI subgroups. We found 
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Table 1 

Demographics and clinical scores for the participants. 





NC 


MCI 


e-MCI 


1-MCI 


AD 


p-value for group difference 






(n = 44) 


(n = 88) 


(n = 62) 


(n = 26) 


(n = 23) 


NC vs MCI 


MCI vs AD 


e-MCI vs I-MCI 


NC vsAD 


Age 


72.7 ± 5.9 


73.3 ± 7.3 


74.0 ± 7.9 


71.8 ± 5.3 


75.8 ± 10.0 


0.59 


0.27 


0.15 


0.17 


Sex 


22 M/22 F 


53 M/35 F 


38 M/24 F 


15M/11 F 


15 M/8 F 










Education 


16.6 ± 2.7 


15.9 ± 2.7 


15.8 ± 2.8 


16.2 ± 2.5 


15.0 ± 3.0 


0.18 


0.22 


0.50 


0.05 


MMSE 


28.9 ± 1.3 


27.8 ± 1.6 


28.0 ± 1.5 


27.5 ± 1.8 


23.0 ± 1.9 


1.01E-4 


4.90E-12 


0.24 


1.57E-14 


ADAS-cog-" 


5.3 ± 2.9 


9.6 ± 4.2 


8.4 ± 3.7 


12.7 ± 4.0 


20.4 ± 7.3 


2.90E-9 


2.04E-6 


1.45E-4 


1.04E-8 


CDR-sob 


0.03 ± 0.1 


1.3 ± 0.7 


1.2 ± 0.6 


1.4 ± 0.9 


4.96 ± 1.4 


138E-27 


2.69E-12 


0.27 


2.68E-14 


ADNI-MEM" 


0.84 ± 0.52 


0.22 ± 0.46 


0.36 ± 0.42 


-0.12 ± 0.38 


-0.74 ± 0.68 


1.01E-8 


4.14E-6 


5.97E-6 


3.64E-10 


ADNI-EF b 


0.79 ± 0.74 


0.16 ± 0.64 


0.14 ± 0.61 


0.19 ± 0.73 


-0.86 ± 0.87 


1.65E-5 


5.16E-5 


0.75 


2.51E-8 



Bold entries signify p < 0.05. 

a ADAS-cog data were available only for a subset of the subjects, with the following numerical breakdown: NC n = 41 , MCI n = 78 (e-MCI = 57. 1-MCI = 21 ), AD n = 20. 
b ADNI-MEM and ADNI-EF composite scores were available only for a subset of the subjects, with the following numerical breakdown: NC n = 41,MCIn = 82 (e-MCI = 58, 1-MCI = 
24), AD n = 20. 



differences in ADAS-cog scores — a measure developed specifically 
for Alzheimer's disease, and the most widely used primary outcome 
measure in clinical trials that test AD drug treatment efficacy (Mohs 
et al., 1983). ADNI-MEM was also significantly different as it is calculated 
using both ADAS-cog and Logical Memory data, the measure used to ini- 
tially classify the two MCI subgroups. 

2.2. MRI and DTI scanning 

All subjects underwent whole-brain MRI scanning on 3 Tesla GE 
Medical Systems scanners at 14 acquisition sites across North America. 
Anatomical Tl -weighted SPGR (spoiled gradient echo) sequences 
(256 x 256 matrix; voxel size = 1.2 x 1.0 x 1.0 mm 3 ; TI = 400 ms; 
TR = 6.98 ms; TE = 2.85 ms; flip angle = 11°), and diffusion- 
weighted images (DWI; 256 x 256 matrix; voxel size: 2.7 x 2.7 x 
2.7 mm 3 ; TR = 9000 ms; scan time = 9 min; more imaging details 
can be found at http://adni.loni.ucla.edu/wp-content/uploads/2010/ 
05/ADNI2_GE_3T_22.0_T2.pdf) were collected. 46 separate images 
were acquired for each DTI scan: 5 T2- weighted images with no diffu- 
sion sensitization (b 0 images) and 41 diffusion-weighted images 
(i> = 1000 s/mm 2 ). This protocol was chosen after conducting a de- 
tailed comparison of several different DTI protocols, to optimize the 
signal-to-noise ratio in a fixed scan time (Jahanshad et al., 2010a; 
Zhan et al., in press). All Tl-weighted MR and DWI images were checked 
visually for quality assurance to exclude scans with excessive motion 
and/or artifacts; all scans were included. 

2.3. Image analysis 

2.3.1. Preprocessing steps 

For each subject, all raw DWI volumes were aligned to the average b 0 
image (DTI volume with no diffusion sensitization) using the FSL 
eddy_correct tool (www.fmrib.ox.ac.uk/fsl) to correct for head motion 
and eddy current distortions. All extra-cerebral tissue was roughly re- 
moved from the Tl-weighted anatomical scans using a number of soft- 
ware packages, primarily ROBEX, a robust automated brain extraction 
program trained on manually "skull-stripped" MRI data (Iglesias et al., 
2011) and FreeSurfer (Fischl et al., 2004). Skull-stripped volumes 
were visually inspected, and the best one selected and sometimes fur- 
ther manually edited. Anatomical scans subsequently underwent inten- 
sity inhomogeneity normalization using the MNI nujcorrect tool (www. 
bic.mni.mcgill.ca/software/). Non-brain tissue was also removed from 
the diffusion-weighted images using the Brain Extraction Tool (BET) 
from FSL (Smith, 2002). To align data from different subjects into the 
same 3D coordinate space, each Tl -weighted anatomical image was lin- 
early aligned to a standard brain template (the downsampled Colin27 
(Holmes et al., 1998): 110 x 110 x 110, with 2 mm isotropic voxels) 
using FSL flirt (Jenkinson et al., 2002) with 6 degrees of freedom (dof) 
to allow translations and rotations in 3D. To correct for echo-planar 



imaging (EPI) induced susceptibility artifacts, which can cause distor- 
tions at tissue-fluid interfaces, skull-stripped b 0 images were linearly 
aligned (FSL flirt, 9 dof) and then elastically registered to their respec- 
tive Tl -weighted structural scans using an inverse-consistent registra- 
tion algorithm with a mutual information cost function (Leow et al., 
2007) as described in (Jahanshad et al., 2010b). The resulting 3D defor- 
mation fields were then applied to the remaining 41 DWI volumes prior 
to estimating diffusion parameters. To account for the linear registration 
of the DWI images to the structural Tl-weighted scan, a corrected gra- 
dient table was calculated. 

2.3.2. DTI maps 

A single diffusion tensor (Basser et al., 1994), or ellipsoid, was 
modeled at each voxel in the brain from the eddy- and EPI-corrected 
DWI scans using FSL dtifit, and scalar anisotropy and diffusivity maps 
were obtained from the resulting diffusion tensor eigenvalues (Ki, \ 2 , 
\ 3 ) which capture the length of the longest, middle, and shortest axes 
of the ellipsoid. The tensor model in DTI has limitations, especially in re- 
gions where fibers cross, but we do not investigate it further here; our 
other papers consider this in more detail (Zhan et al., 2008; Leow 
et al., 2009; Zhan et al., 2009; Zhan et al., 2010; Zhan et al., 2011 ; Zhan 
et al., 2012b; Zhan et al., in press). Fractional anisotropy (FA), a measure 
of the degree of diffusion anisotropy, was calculated from the standard 
formula: 

pA= fi y/fa -<\>) 2 + (\ 2 -<\.>) 2 + (\ 3 -<k>) J 
V2 ^\ 2 + \ 2 + \ 2 

<x > _ + X 2 + ^3 



Table 2 

Index of 43 ROls from the WM tract atlas (Mori et al., 2008) followed by their 
abbreviations. 



Genu of corpus callosum 


GCC 


Posterior thalamic radiation 


PTRL.R 


Body of corpus callosum 


BCC 


Sagittal stratum 


SSL.R 


Splenium of corpus callosum 


see 


External capsule 


EC L.R 


Full corpus callosum 


cc 


Cingulum (cingulate gyrus) 


CGC L,R 


Corticospinal tract 


CSTL.R 


Cingulum (hippocampus) 


CGHL.R 


Cerebral peduncle 


CPL.R 


Fornix (crus)/Stria terminalis 


FX/ST L,R 


Anterior limb of internal 


ALIC L,R 


Superior longitudinal 


SLFL.R 


capsule 




fasciculus 




Posterior limb of internal 


PLIC L,R 


Superior fronto-occipital 


SFOL.R 


capsule 




fasciculus 




Retrolenticular part of internal 


RUC L.R 


Inferior fronto-occipital 


IFO L,R 


capsule 




fasciculus 




Anterior corona radiata 


ACRL.R 


Uncinate fasciculus 


UNC L,R 


Superior corona radiata 


SCRL.R 


Tapetum 


TAP L,R 


Posterior corona radiata 


PCRL.R 


All ROls 


TOTAL 
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- where <\> is the mean diffusivity (MD), or average rate of diffusion in 
all directions. Axial diffusivity was defined as the primary (largest) 
eigenvalue (AxD = and captures the longitudinal diffusivity, or 
the diffusivity parallel to axonal fibers (assuming of course that the 
principal eigenvector is indeed following the dominant fiber direction, 
which may be unclear in regions with extensive fiber crossing). Radial 
diffusivity (RD), which captures the average diffusivity perpendicular 
to axonal fibers, was calculated as the average of the two smaller eigen- 
values: 



2.3.3. White matter tract atlas ROI summary measures 

We linearly, then elastically registered (Leow et al., 2007) the FA 
image from the Johns Hopkins University (JHU) DTI atlas (Mori et al., 
2008) to each subject's distortion corrected FA image. We then applied 
that deformation to the stereotaxic JHU "Eve" atlas WM labels (http:// 
cmrm.med.jhmi.edu/cmrm/atlas/human_data/file/AtlasExplanation2. 
htm), using nearest neighbor interpolation to avoid intermixing of la- 
bels. This is not the atlas that had the problem pointed out by Rohlfing 
(2013). We then superimposed the atlas ROls into the same coordinate 
space as our results. We removed 10 ROls from the analyses (including 
the left and right middle cerebellar peduncle, pontine crossing tract, 
medial lemniscus, inferior and superior peduncles) as they often fell 
partially or completely out of the field-of-view (FOV) of the images. 
We also excluded the body of the fornix as it is small and prone to 
misregistration and partial voluming. AD researchers are specifically in- 
terested in the fornix as it is the primary posterior pathway coming out 
of the back of the hippocampus, a key target of pathology. While we in- 
cluded the crus of the fornix/stria terminalis, the body is just too small to 
be resolved well on DTI at this resolution. In addition to the JHU labels, 4 
more ROls were evaluated: bilateral genu, body, and splenium of the 
corpus callosum (as opposed to just the lateralized measures), as well 
as the entire corpus callosum, and a large "TOTAL" WM ROI made up 
of all the other ROls, to obtain total summary measures of these regions. 
We were then able to calculate the average FA, MD, RD and AxD, within 
the boundaries of each of the 43 ROls for each subject (Table 2). 2 

2.3.4. TBSS tract atlas ROI summary measures 

Tract-based spatial statistics (TBSS) (Smith et al., 2006), provided in 
the FSL software package (http://www.fmrib.ox.ac.uk/fsl/), was also 
performed according to protocols outlined by the ENIGMA-DTI group: 
http://enigma.loni.ucla.edu/ongoing/dti-working-group/. All subjects' 
corrected FA maps were linearly, then elastically registered (Leow 
et al., 2007) to the ENIGMA-DTI template in ICBM space. The resulting 
3D deformation fields were then applied to the three diffusivity maps. 
All subjects' spatially normalized FA, MD, RD and AxD data were 
projected onto the skeletonized ENIGMA-DTI template. Mean anisotro- 
py and diffusivity measures were calculated along the skeleton in the 
same 43 ROls (Table 2). This type of analysis has been used previously 
in both genetic studies and studies of disease to home in on associated 
WM tracts (Kochunov et al., 201 1 ; Jahanshad et al., 2013). 

2.3.5. Template creation and spatial normalization 

A study-specific minimal deformation template (MDT) (Gutman 
et al., 2012) was created using 29 cognitively healthy elderly control 
(NC) spatially aligned FA maps. An MDT deviates, on average, the least 
(in some metric) from the anatomy of the subjects, and can often im- 
prove registration accuracy and statistical power (Gutman et al., 2012; 
Lepore et al., 2007). The MDT was generated by creating an initial affine 
mean template from all 29 subjects, then registering all the aligned 



2 Except for the total WM ROI, all these ROI measures for all subjects are available for 
download at https://ida.loni.ucla.edu/. 



individual scans to that mean using a fluid registration (Leow et al., 
2007) while regularizing the Jacobians (Yanovsky et al., 2007). A new 
mean was created from the registered scans; this process was iterated 
several times. Each subject's initial FA map was elastically registered to 
the final MDT and the resulting deformation fields were applied to the 
3 diffusivity maps to align them to the same coordinate space. To ensure 
white matter alignment across subjects, registered FA maps were 
thresholded at FA > 0.2 to include only highly anisotropic anatomy and 
the thresholded maps were elastically registered to the thresholded 
MDT (FA > 0.2). Again, the resulting deformation fields were applied to 
all previously registered DTI maps. We also used the tissue-specific, 
smoothing-compensated method {T-SPOON) proposed by Lee et al. 
(2009) to improve tissue specificity and reduce confounds caused by 
morphometric differences that are not fully corrected by the elastic 
registration. 

2.4. Statistical analysis 

We ran voxel-wise multiple linear regressions, covarying for age and 
sex, to test for statistical effects of AD and MCI diagnosis relative to the 
NC group, on measures of white matter integrity in FA, MD, RD, and AxD 
maps. We also tested for associations between these DTI measures and 
MMSE, CDR-sob, and ADAS-cog scores, controlling for age and sex, 
across the entire study sample (i.e., including all AD, MCI and NC sub- 
jects), and also within each subgroup. We also ran all of the regressions 
using a random-effects regression model, grouping the data by acquisi- 
tion site. To limit statistical testing to white matter, where the power is 
greater to detect differences, statistics were run only on voxels within 
the boundaries of the MDT thesholded at FA > 0.2. Prior studies have 
thresholded FA values between 0.2 and 0.3 to successfully exclude 
gray matter or CSF (Wakana et al., 2005; Smith et al., 2006). As we 
were studying a clinical population, we chose the more conservative 
(lower) limit of the recommended FA threshold. We also only ran statis- 
tics on voxels of the thresholded MDT present in all subject scans, as 
some scans had a slightly cropped FOV. As such, we did not consider 
the inferior parts of the cerebellum and brain stem. 

We ran random-effects regressions on the average anisotropy and 
diffusivity measures within the 43 full ROls, again covarying for age 
and sex, testing for statistical effects of diagnosis (NC vs. MCI, or AD), 
3 global clinical test scores, and post-hoc tests on ADNI-MEM and 
ADNI-EF. We further tested and compared TBSS ROI measures. 

Computing multiple association tests for each ROI, or thousands of 
tests on a voxel-wise level can introduce a high false positive error 
rate. To account for these errors, we used the standard false discovery 
rate (FDR) method to control the false positive rate of each map at 
q = 0.05 (Benjamini and Hochberg, 1995). All statistical maps shown 
in this paper (Figs. 1-5) were thresholded at the FDR critical p-value. 

To visualize and compare effect sizes in the anisotropy and diffusiv- 
ity maps, we computed and graphed the cumulative distribution func- 
tions (CDF) of the p-values obtained from the voxel-wise and ROI 
random-effects regressions. The ordered set of observed p-values was 
plotted against the expected null distribution. If there are no group dif- 
ferences (a null distribution), then the plot would fall approximately 
along the y = x diagonal line. However, if the CDF initially rises at 
a rate steeper than 20 times the null CDF (y = 20x), then the corre- 
sponding maps are FDR significant at q = 0.05. A greater slope reflects 
a larger effect size; the FDR critical p-value (at q = 0.05) is the point at 
which the curve intersects they = 20x line. 

DTI measures calculated from a single-tensor diffusion model have 
limitations in regions with extensive fiber crossing. For example, in AD 
we might expect FA deficits throughout the brain, but FA may appear 
to be artificially increased where crossing fibers deteriorate (Douaud 
et al., 201 1 ). To make sure that these results were not false positives 
we separated the full p-value map (irrespective of the significance 
level of the voxels ) into two maps : ( 1 ) voxels that showed a negative as- 
sociation, and (2) positively associated voxels. We then independently 
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corrected each p-map for multiple statistical comparisons using FDR 
(Benjamini and Hochberg, 1995). 

As an alternative way to visualize white matter integrity differences 
between patients and controls, we also created maps showing the "per- 
cent difference" in mean DTI measures between AD patients and NC, 
within the FDR significant regions. 



There is a mild bias in doing this, in that it selects voxels that show 
effects, and then computes the magnitude of the effects. Even so, it 
was used to simply describe how much difference was seen in the re- 
gions shown, bearing in mind that the same difference would not be 
found in the rest of the brain. 

3. Results 

3.1. Voxel-based analyses 

3.1.1. White matter integrity differences between diagnostic groups 

3.1. 1.1. AD vs controls. As measured by FA, AD patients showed pervasive 
deficits in fiber integrity (lower FA) compared to healthy controls, 
throughout the WM across the entire brain, when using a general linear 
regression model to adjust for age and sex (critical p < 0.016; 31.9% of 
voxels survived the FDR threshold (Benjamini and Hochberg, 1995); 
minimum p-value, 1.65 x 10~ 8 ). However, when we used a random- 
effects regression model at every voxel to take into account acquisition 
site differences, we found more significant and even more widespread 
results (critical p < 0.020; 40.1% of voxels survived the FDR threshold 
(Benjamini and Hochberg, 1995); minimum p-value, 6.66 x 10~ 9 ; 
Fig. la). We therefore proceeded with this model - fitting a site effect - 
for the remainder of the analyses. A few significant voxels exhibited as- 
sociations in a direction opposite to the great majority of the brain map, 
and contrary to what would traditionally be accepted as showing impair- 
ment. For example, we found regions where the AD group had a higher 
mean FA (Fig. le; blue regions). These regions were largely found at the 
junction of the corpus callosum commissural fibers and the corona 
radiata. Such regions are notorious for fiber crossings, which may artifi- 
cially reduce FA, if FA is estimated using the single-tensor diffusion 
model (Oishi et al, 2011). This pattern has been reported in other AD 
studies, and may reflect a selective sparing or selective degeneration of 
one of the pathways in a region with crossing fibers (Douaud et al., 
2011). When we separated (and independently FDR-corrected) the 
map of voxels that showed a negative association with AD and the 
map of the positively associated voxels, the diffuse negative associations 
were still significant, but the previously significant positive associations 
were not. 

This pervasive pattern of significance was largely replicated for 
the three other diffusivity measures. As predicted, we found higher 
mean RD, AxD and MD values in AD patients relative to controls 
(Fig. 1 b-d) and the effects were larger and even more diffuse than for 
FA. CDF plots confirm that FA was not the most sensitive measure for 
differentiating groups. Comparisons of NC to AD patients (Fig. 2a) re- 
vealed that RD was the most sensitive measure followed closely by MD. 

3.1.1.2. Early, late and all MCI vs controls. No significant difference was 
detected between e-MCI and 1-MC1 groups for any of the DTI measures, 
so we first assessed the MCI group as a whole, followed by these sub- 
groups. We found no significant differences between NC and MCI as a 
whole (n = 88), or NC and the e-MCI (n = 62) subgroup, for any of 
the anisotropy or diffusivity measures. However, for NC compared to 
1-MCI subjects (n = 26), we did find significantly higher diffusivity 
measures in the left hippocampal part of the cingulum (MD: FDR critical 
p < 3.05 x 10~ 5 ; RD: FDR critical p < 2.72 x 10~ 5 ; AxD: FDR critical 



p < 2.46 x 10 -5 ; Fig. 3). While these measures are significant on their 
own, the N was small. These contrasts will benefit from further analysis, 
as the ADN1 dataset grows. 

3.3.2. Correlation with MMSE, CDR-sob, and ADAS-cog neuropsychological 
scores 

3. 1.2.1. Full cohort. We first assessed anisotropy and diffusivity map asso- 
ciations with widely used clinical cognitive ratings in the entire study 
population, including AD patients, MCI, and NC subjects (Fig. 4). Fig. 4a 
shows WM regions where anisotropy and diffusivity differences corre- 
lated with the MMSE scores in the entire study population (n = 155). 
MMSE was significantly positively associated with FA (FDR critical 
p < 0.018), and negatively associated with diffusivity (MD: FDR critical 
p < 0.037; RD: FDR critical p < 0.040; AxD: FDR critical p < 0.028). 
That is, lower FA and higher diffusivity, which typically indicate greater 
WM deficits, were associated with lower MMSE scores, indicative of 
greater impairment. Fig. 4b shows WM regions where FA and diffusivity 
differences correlated with the CDR-sob scores (n = 155). As expected, 
FA was negatively associated with CDR-sob (FDR critical p < 0.027), and 
diffusivity was positively associated (MD: FDR critical p < 0.039; RD: 
FDR critical p < 0.042; AxD: FDR critical p < 0.029). Fig. 4c shows WM 
regions where FA was negatively associated with ADAS-cog (n = 139; 
FDR critical p < 0.016), and diffusivity positively associated (MD: FDR 
critical p < 0.040; RD: FDR critical p< 0.041; AxD: FDR critical 
p < 0.029). Lower FA and higher diffusivity were associated with higher 
(worse) ADAS-cog and CDR-sob scores. CDF plots of MMSE, ADAS-cog, 
and CDR-sob voxel-wise associations again reveal that FA associations, 
while FDR significant, show the smallest effect size, while MD and RD 
have the largest effect sizes (Fig. 2d). 

As with diagnostic group differences, in all clinical score analyses we 
found small significant regions with associations in a direction opposite 
to what would traditionally be accepted as showing impairment. Again, 
when we corrected for multiple comparisons on voxel-wise p-maps of 
each association direction separately, the associations in the unexpected 
direction were no longer significant and only the intuitive direction 
remained. 

3.1.22. Early, late and all MCI. We further assessed associations between 
cognitive scores and DTI measures in just the MCI group, as well as 
e-MCI and 1-MCI subgroups. We found a significant positive association 
between all three diffusivity measures and CDR-sob within the MCI 
group (MD: FDR critical p < 0.006; RD: FDR critical p < 0.006; AxD: 
FDR critical p < 0.002; Fig. 5) and a small positive association between 
ADAS-cog and MD in the middle or lateral occipital WM in 1-MCI (FDR 
critical p < 5.35 x 10 -6 ). However, there was no detectable association 
with anisotropy or diffusivity measures and MMSE scores. 

3.2. ROI analyses 

32.1. ROI differences between diagnostic groups 

3.2.1.1. AD vs controls. We sorted the JHU atlas ROI results by p-value to 
assess which ROIs showed the greatest differences between groups. In 
doing so, we have to bear in mind that this is not a ranking of the disease 
effects on the brain by severity; results may also be affected by how ac- 
curately each region can be measured, which in turn depends on the 
size and the homogeneity of the region. With these caveats in mind, 
Tables 3-4 list the 'top 10' most significant ROIs for each regression 
analysis (here and in what follows, we use the term 'most significant' 
to mean lowest p-values in the tests of group differences in mean 
values; we acknowledge that this term is not always preferred by stat- 
isticians as the p-values from different tests cannot necessarily be com- 
pared). When comparing NC to AD patients (Table 3a), the splenium of 
the corpus callosum, left fornix [crus)/stria terminalis, left hippocampal 
part of the cingulum, and total ROI map, were in the top 10 most 
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Fig. 1. (Top panels) Statistical maps show — log 10 p-values within regions where AD patients exhibit significantly higher (a) MD (FDR critical p < 0.038), (b) RD (FDR critical p < 0.041). and 
(c) AxD (FDR critical p < 0.027) than healthy controls (NC). The cingulum, temporal lobe (including the hippocampal part of the cingulum) and splenium of the corpus callosum show the 
"most significant" differences, i.e. greatest effect sizes in the maps, (d) Significant but slightly less profuse FA differences between groups were also found (FDR critical p < 0.020). (e) Beta- 
values (non-normalized slope of the regression) within FDR significant regions largely revealed regions with lower FA in AD, but there are small regions, notorious for crossing fibers, 
where AD patients have higher mean FA than controls (dark blue), contrary to what would traditionally be accepted as showing impairment. (Bottom panel) Maps show the percent dif- 
ference in FA and diffusivity measures between AD patients and NC within FDR significant regions. Again, the corpus callosum, temporal lobe, cingulum ( including the hippocampal part of 
the cingulum), and the posterior thalamic radiations show the greatest degree of difference (up to -33%); we note that the significance values and effect sizes in these regions depend on 
the mean group differences and its standard error, which is more completely reflected in top panels. In both figures, MD and RD differences are greater and more widespread. (For inter- 
pretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 
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Fig. 2. Effect sizes for FA, MD, RD and AxD measures in AD, MCI and NC group comparisons and clinical cognitive test score associations reveal that FA is consistently the least sensitive 
measure. This type of plot is used to compare effect sizes in statistical maps based on different diffusivity measures. We used it to help decide which DTI measures best distinguish the 
diagnostic groups. Here we show cumulative distribution function (CDF) plots of the distribution of the p-values obtained from (a,d) voxel-wise and (b-c, e) ROI linear regressions, 
which are subjected to multiple comparisons correction using standard FDR (Benjamini and Hochberg, 1995). (a) VBA comparisons of NC to AD patients revealed that RD was the 
most sensitive measure (denoted by the higher critical p-values controlling the FDR, i.e., the highest non-zero x-coordinate where the CDF crosses the y = 20x line) followed closely 
by MD. (b) ROI comparisons of NC vs AD confirm that FA was the least sensitive measure, (c) ROI comparisons of NC vs MCI reveal that AxD measures had the largest effect size, followed 
by MD (VBA comparisons were not significant), (d) Voxel-wise and (e) ROI linear regressions on associations between cognitive scores (MMSE, CDR-sob, and ADAS-cog) in the entire 
cohort and anisotropy and diffusivity measures confirm that FA associations, while FDR significant, show the smallest effect size, while MD and RD have the largest effect sizes. 



significant ROIs, across all anisotropy and diffusivity measures. While 
the superior corona radiata and left sagittal stratum (which includes in- 
ferior longitudinal fasciculus and inferior fronto-occipital fasciculus) 
were in the top 10, across all diffusivity measures. 

3.2. 1.2. Early, late and all MCI. As there were no significant differences be- 
tween 1-MC1 and e-MCI subjects across any of the measures, we com- 
pared the entire MCI group to NC subjects (Table 3b). While no ROI 
showed significant FA differences, the splenium of the corpus callosum, 
left tapetum, left hippocampal part of the cingulum, and left fornix 
(cms)/stria terminalis were in the top 1 0, among all diffusivity measures. 
We further assessed differences between e-MCl and NC and 1-MC1 and 
NC. The left hippocampal part of the cingulum was significant across 
all three diffusivity measures in both comparisons (Table 3c-d). FA 
showed no statistically significant differences. 

CDF plots of the ROI p-values for NC vs AD (Fig. 2b) confirm VBA 
findings that FA was the least sensitive measure, and that RD has a 
slightly higher effect size than other measures. However, CDF plots of 
the ROI p-values, for NC vs MCI, reveal that AxD measures had the larg- 
est effect size, followed by MD (Fig. 2c). Table 3d also reveals that AxD 
has a larger effect size when comparing 1-MCI and NC. 



3.2.2. ROI correlations with clinical global neuropsychological scores 

The strongest associations with MMSE score (Table 4a) across all 
four measures were found in the left hippocampal part of the cingulum, 
the left fornix (crus)/stria terminalis, and the total ROI map. All three dif- 
fusivity measures were additionally associated with MMSE in the left 
cingulum, and the bilateral sagittal stratum (this includes the inferior 
longitudinal fasciculus and inferior fronto-occipital fasciculus). 

The strongest associations with CDR-sob score (Table 4b) across all 
four measures were found in the entire corpus callosum, the left hippo- 
campal part of the cingulum, and the total ROI map. All three diffusivity 
measures additionally associated with CDR-sob in the left cingulum, and 
bilaterally in the sagittal stratum and inferior fronto-occipital fasciculus. 

The strongest associations with the ADAS score (Table 4c) across all 
four measures were found in the left hippocampal part of the cingulum, 
and in all three diffusivity measures in the left cingulum, bilaterally in 
the sagittal stratum and inferior fronto-occipital fasciculus, and in the 
total ROI map. 

As with voxel-wise measures, we further assessed the association 
between cognitive scores and ROI DTI measures in just the MCI group 
as a whole and also in the e-MCl and 1-MCI groups separately 
(Table 5). As would be predicted, we found significant negative MMSE 
associations with diffusivity in both the MCI group and e-MCI subgroup 
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with ADNI-MEM bilaterally in the cingulum and hippocampal cingulum, 
bilateral inferior fronto-occipital fasciculus and sagittal stratum, and left 
uncinate and external capsule. 

These measures were not linked to DTI measures in the MCI group as 
a whole, nor the e-MCI and 1-MCI groups separately. 

4. Discussion 




5.5 6.0 

|-log 10 (p-value) | 

Fig. 3. Statistical maps show the — log 10 p-values within regions where 1-MCI subjects ex- 
hibit significantly higher AxD (FDR critical p < 2.46 x 10~ 5 ) than healthy controls (NC). 
The left hippocampal part of the cingulum was also significant in MD and RD maps. This 
region is small, but it passes the conventional FDR correction for multiple comparisons, 
and is in a region implicated in early changes in AD (medial temporal regions). With a 
greater sample size, this effect may be detectable in the other hemisphere as well. 



in the left hippocampal part of the cingulum. We found significant pos- 
itive associations between MD and AxD and CDR-sob scores in the MCI 
group and no associations with ADAS-cog. 

Again, CDF plots of the ROI p-values for cognitive score associations 
(Fig. 2b) confirm that FA was the least sensitive measure. 

32.3. ROI analyses summary 

ROIs that were consistently most significant across all three diffusiv- 
ity measures in more than one analysis (Fig. 6) included the splenium of 
the CC, the left cingulum, particularly the hippocampal part, the left 
fornix (crus)/stria terminalis which projects to the dorsal region of the 
hippocampus, bilateral temporal lobe sagittal stratum, and the average 
measure across all ROIs. They also included the bilateral inferior 
fronto-occipital fasciculus in clinical score analyses. Of these ROIs, the 
left fornix (crus)/stria terminalis and hippocampal cingulum, and total 
ROI map were most significant across all anisotropy and diffusivity mea- 
sures in at least two analyses. However, the left hippocampal part of the 
cingulum was in the top 10 ROIs in by far the most analyses. 

For complete tables of TBSS ROI diagnostic group and cognitive score 
analyses results, please see Inline Supplementary Tables S2-S4. To sum- 
marize, TBSS ROIs that were consistently most significant (i.e., giving 
greatest effect sizes) across all four anisotropy and diffusivity measures 
in more than one analysis included the splenium of the CC, bilateral hip- 
pocampal part of the cingulum, the left sagittal stratum and uncinate, 
and the average measure across all ROIs (Inline Supplementary 
Fig. SI ). While full ROI analyses additionally revealed bilateral inferior 
fronto-occipital fasciculus, left fornix (crus)/stria terminalis, left cingu- 
lum, and right sagittal stratum, as top ROIs and TBSS did not, TBSS 
ROIs additionally revealed the left uncinate and right hippocampal 
part of the cingulum. For a comparison of TBSS and full ROI analyses 
p-value CDF plots see Inline Supplementary Fig. S2. 

Inline Supplementary Tables S2-S4, Figs. SI and S2 can be found on- 
line at http://dx.doi.Org/10.1016/j.nicl.2013.07.006. 

3.2.4. Post-hoc ROI correlations with memory and executive function scores 
The strongest associations between full ROI measures and ADNI-EF 
(Table 6a) in the entire study cohort (n = 143) across all four DTI mea- 
sures were found in the bilateral sagittal stratum and the total ROI map. 
All three diffusivity measures were additionally associated with ADNI- 
EF in the bilateral cingulum, left hippocampal cingulum, and left superi- 
or longitudinal fasciculus. 

There were no associations between FA and ADNI-MEM (Table 6b) 
in the entire cohort, but all three diffusivity measures were associated 



In this ADNI study, we found that DTI indicators of white matter im- 
pairment have the potential to emerge as useful clinical tools for differ- 
entiating diagnostic groups in studies of AD. We had three main 
findings: ( 1 ) in ROI and VBA analyses, we found widespread anisotropy 
and diffusivity disruptions, often in tracts that pass through the tempo- 
ral lobe and posterior brain regions ( especially the left hippocampal cin- 
gulum), in elderly AD and MCI patients; (2) these disruptions were also 
associated with neuropsychological and cognitive deficits; and (3) dif- 
fusivity (MD, RD, AxD) measures were more sensitive for detecting dif- 
ferences than FA measures, and could detect more subtle MCI 
differences, where FA could not. 

Alzheimer's disease (AD) is characterized by neuronal loss and wide- 
spread gray matter atrophy, but it is also marked by a disturbance in the 
brain's WM pathways, perhaps secondary to the effect of cortical neuro- 
nal loss. Changes in white matter neuropathology include partial loss of 
axons, myelin sheaths, and oligodendroglial cells (Brun and Englund, 
1986; Sjobeck et al., 2005). AD patients have been shown to have signif- 
icantly more WM hyperintensities (WMH) than controls. WMH are sig- 
nificantly related to cortical atrophy (Capizzano et al., 2004), entorhinal 
cortex (Guzman et al., in press) and hippocampal atrophy (de Leeuw 
et al., 2004) in AD patients. Significant WM atrophy has also been re- 
ported (Hua et al., 2008; Hua et al., 2010; Migliaccio et al., 2012). 

There is growing diffusion imaging evidence of AD related WM 
changes as well. Most studies report lower FA and higher MD in all 
lobes of the brain in both MCI and AD when compared to cognitively 
healthy controls, with emphasis on medial temporal lobe structures 
(Kantarci et al., 2001; Bozzali et al., 2002; Takahashi et al., 2002; 
Fellgiebel et al., 2004; Medina et al., 2006; Rose et al., 2006; Xie et al., 
2006; Zhang et al., 2007; Kavcic et al., 2008; Stebbins and Murphy, 
2009), consistent with regions showing earliest pathological changes 
(Braak and Break, 1991; Braak and Braak, 1995; Thompson et al., 2007) 
and our results. Many DTI methods have been used to assess WM differ- 
ences in AD including TBSS (Smith et al., 2006) and tract or connectivity 
analyses (Daianu et al., 2012, 2013a,b). Here we use the most common 
methods, VBA and ROI analyses, with ROIs based on a stereotaxic WM 
ROI atlas (Mori et al., 2008), which has been shown to accurately 
parcellate anatomical regions in AD patients with severe atrophy 
(Oishi et al., 2009). Using predetermined, template-based ROIs can 
limit findings, and registration across subjects may vary due to morpho- 
logical differences. However, using a common atlas helps to make quan- 
tification efficient and easier to standardize across sites (Jahanshad et al., 
2013). Given that VBA findings are diffuse, ROIs are arguably helpful to 
begin to rank the most sensitive regions for detecting group differences. 

When comparing AD patients to NC, the largest VBA differences 
were found in the corpus callosum (CC), temporal lobe, cingulum and 
hippocampal cingulum, and regions near the posterior thalamic radia- 
tions. ROI analyses revealed similar diffuse patterns with the CC 
splenium, the left fornix {crus)/stria terminalis, and average measure 
across all ROIs in the top 10 most significant (lowest p-values) ROIs 
across all anisotropy and diffusivity measures. However, we found 
that the most significant ROIs were not necessarily the most coherent 
tracts. Large clinical DTI studies with limited spatial and directional res- 
olution are often better powered to find deficits in regions where the FA 
and fiber coherence is highest, such as the splenium of the CC. In a sup- 
plementary test, we assessed whether the ROI average FA values in the 
NC group was correlated with the effect size of each ROI (which is relat- 
ed to the p-values obtained) in the NC vs AD comparison (see supple- 
mentary text). Inline Supplementary Fig. S3 shows that the ROI 
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Fig. 4. Statistical maps show that clinical scores on the (a) MMSE, (b) CDR-sob, and (c) ADAS-cog are related to detectable differences in fractional anisotropy, and mean, radial and axial 
diffusivity in the entire cohort. These maps show the — log 10 p-values within regions that significantly correlate with cognitive scores. RD and MD show the most widespread and robust 
associations. The hippocampal part of the cingulum and surrounding temporal lobes show the "most significant" associations (i.e., lowest p-values and highest voxel-wise effect sizes) 
between cognitive performance and white matter integrity. 



average FA values in the NC group were not correlated with the effect 
size (defined by the Z-score corresponding to the p-values) of each ROI. 

Only diffusivity measures were sensitive enough to detect more sub- 
tle differences between NC and MCI groups. ROls including the 
splenium of the CC, left tapetum, left hippocampal part of the cingulum 
and left fornix (crus)/stria terminalis consistently differed in all diffusiv- 
ity measures but not FA. Only the left hippocampal part of the cingulum 
was significant across all three diffusivity measures in both comparisons 
of e-MCl and 1-MCI to controls. Similarly, only the left hippocampal part 
of the cingulum displayed significant increased diffusivity in 1-MCI in 
the VBA analyses. 

Inline Supplementary Fig. S3 can be found online at http://dx.doi. 
org/1 0.1 01 6/j.nicl.201 3.07.006. 

The lack of larger differences between the MCI and normal control 
groups is surprising in light of other studies showing that these groups 
often differ on other neuroimaging measures. This may be due the small 



sample size and also heterogeneity in the MCI cohort. The MCI group re- 
ferred to in papers using ADNI1 data is now called late MCI (1-MCI) in 
ADNI2. In ADNI2, the majority of the MCI group includes the enrollment 
of a new cohort called early MCI (e-MCI), with milder episodic memory 
impairment than the 1-MCI group. Among the MCI subjects with DTI 
data available, a smaller percent will be 1-MCI than in phase one of 
ADNI, perhaps contributing to the apparent discrepancy with other 
ADNI studies/results. This difference in the ADNI2 versus ADNI1 and 
other MCI cohorts should be considered in expectations of what the ef- 
fect sizes should be. 

Strongest ROI associations with cognitive scores in the full cohort 
were consistently found in the left hippocampal part of the cingulum 
across all anisotropy and diffusivity measures. The total ROI map was 
also one of the top 10 lowest p-values in almost all FA and diffusivity 
analyses, while left cingulum and bilateral sagittal stratum were in the 
top 10 in all diffusivity analyses. The bilateral inferior fronto-occipital 
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Fig. 5. Statistical maps show where CDR-sob scores are significantly positively associated with diffusivity measures in the MCI group, considered on its own (n = 88). Higher diffusivity, 
indicative of greater WM deficits, is associated with greater CDR-sob scores indicative of greater impairment. We find the "most significant" association ( i.e., strongest voxel-wise effect) in 
the splenium of the corpus callosum and the posterior cingulum. 



fasciculus was in the top 10 in almost all full cohort diffusivity analyses. 
When it came to picking up subtle differences within the MCI group, 
CDR-sob scores were most widely correlated with diffusivity in both 
VBA and ROl analyses. Prior studies also suggest that CDR-sob is a more 
sensitive clinical assessment than ADAS-cog and MMSE, and it may better 
relate to measures of atrophy on anatomical MR] (Hua et al„ 2009). 

Post-hoc analyses were conducted to assess if there was any associ- 
ation between DTI summary measures and executive function (ADNI- 
EF) and memory (ADNI-MEM) scores. These analyses revealed some 
of the strongest associations in temporal lobe tracts that subserve a va- 
riety of functions known to deteriorate during aging, including the bilat- 
eral sagittal stratum, that carries parts of the inferior fronto-occipital 
fasciculus and inferior longitudinal fasciculus, the hippocampal cingu- 
lum as well as the cingulum. ADNI-EF associations were additionally 
found in the total ROl and superior longitudinal fasciculus, while 
ADNI-MEM was associated with DTI measures in the external capsule, 
inferior fronto-occipital fasciculus and uncinate fasciculus. These re- 
gional differences could shed light on the network of brain regions, 
connected via white matter fiber bundles, involved in each task. Prior 
DTI studies have linked similar regions to performance in executive 
function and memory (Kantarci et al., 2011; Sasson et al., 2013). 

As patterns of differences were often diffuse across the VBA WM 
map, the total WM ROl was sometimes one of the 10 lowest p-values, 
but it was never the absolute best. Therefore assessing individual ROIs 
may offer slightly more power to detect changes than a global summary 
measure (Jahanshad et al., 2013). In general, the most significant ROIs 
were found mainly in the left hemisphere (for example, the left hippo- 
campal part of the cingulum was in the top 10 ROIs in by far the most 



analyses), which has also been found in some prior MRI and DTI studies 
(Fox et al., 1996; Thompson et al„ 2001 ; Scahill et al., 2002; Thompson 
et al., 2003; Muller et al., 2005). A supplementary analysis revealed no 
evidence of greater variability in ROl measures in the right hemisphere 
that may be causing the better performance of left hemisphere ROIs (see 
supplementary text). TBSS ROl results did not differ dramatically from 
the full average ROl results and left hippocampal cingulum was again 
'top 10' in the most analyses; however, it was closely followed by the 
right hippocampal cingulum. 

Aside from the total ROl and bilateral IFO, the ROIs that were consis- 
tently significant in at least two analyses (Fig. 6) corroborate a pattern of 
degeneration in the temporal lobe and posterior temporo-parietal cir- 
cuitry found in many other DTI studies of MCI and AD (Head et al., 
2004; Stahl et al., 2007; Chua et al., 2008; Stebbins and Murphy, 
2009) ; these include limbic tracts in the parahippocampal white matter, 
posterior cingulum, fornix, and splenium of the CC, and have been 
linked to lower cognitive scores (Rose et al., 2000; Bozzali et al., 2002; 
Takahashi et al., 2002; Yoshiura et al., 2002; Fellgiebel et al., 2004; 
Fellgiebel et al., 2005; Duan et al., 2006; Medina et al., 2006; Rose 
et al., 2006; Zhang et al., 2007; Medina and Gaviria, 2008; Mielke 
et al., 2009). DTI studies have found lower anisotropy in the white mat- 
ter pathway of the cingulate gyrus, particularly the posterior cingulum, 
which connects to the entorhinal cortex and plays a role in the choliner- 
gic system, known to be impaired in AD (Perry, 1980; Selden et al., 
1998; Takahashi et al., 2002; Medina et al., 2006; Zhang et al., 2007). 
In fact, it has been highly implicated in studies comparing MCI to con- 
trols, especially on the left as in our study (Fellgiebel et al., 2005; 
Medina et al., 2006; Rose et al., 2006; Zhang et al., 2007; Chua et al., 
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Table 3 

Atlas ROI average anisotropy and diffusivity. While many regions were statistically significant (Benjamini and Hochberg, 1995), here we highlight regions with the "top 10" significant p- 
values (greatest effect sizes) when comparing mean values between diagnostic groups ( in the left column we show the diagnostic groups being compared ). We also note the total number 
of significant ROIs (after FDR correction) out of the 43 tested and the critical FDR p-value. 
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6.19E-11 


SS L 


6.29E— 1 1 


NC 


sec 


1.40E-6 


CGHL 


1.23E-9 


CGHL 


1.69E-09 


CGHL 


1.68E-9 


vs 


TAP L 


2.46E-6 


TOTAL 


4.63E-9 


TOTAL 


6.89E-09 


TOTAL 


7.47E-9 


AD 


PTRL 


2.07E-5 


SCRL 


1.11E-8 


SCRL 


6.22E-08 


SCRL 


1.74E-8 




CC 


4.41 E- 5 


sec 


6.54E-8 


PTRL 


7.64E-08 


IFOL 


3.04E-8 




SFOR 


5.86E-5 


Fx/STL 


1.07E-7 


Fx/STL 


7.70E-08 


SSR 


9.02E-8 




TAPR 


8.79E-5 


PTRL 


1.10E-7 


sec 


8.05E-08 


Fx/ST L 


2.54E-7 




ACRL 


1.35E-4 


SSR 


1.24E-7 


TAP L 


1.25E-07 


ECL 


3.21 E- 7 




TOTAL 


1.62E-4 


TAP L 


1.43E-7 


ACRL 


1.40E-07 


sec 


334E-7 




CGHL 


1.84E-4 


1FOL 


2.77E-7 


SCRR 


2.82E-07 


UNCL 


3.46E-7 
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# sig/43 
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d) - - CGHL 4.71E-5 CGHL 9.53E-5 CGHL 3.01E-5 

NC - - - - - - SSL 1.88E-4 

vs - - - - - - SSR 1.47E-3 
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#sig/43 0 l(p<4.71E-5) l(p<9.53E-5) 5 (p < 5.678E-3) 



2008). Cognitive scores like MMSE have also been linked to DTI diffusiv- 
ity values in the posterior cingulate gyrus in AD (Yoshiura et al., 2002; 
Fellgiebel et al., 2005). Most AD studies also find WM deficits in the cor- 
pus callosum, but it is not clear whether the genu is more affected (Xie 
et al., 2006), in line with the hypothesis that later maturing regions are 
the first affected, or the splenium, corroborating a pattern of degenera- 
tion in the posterior temporal parietal circuitry (Takahashi et al., 2002; 
Medina et al., 2006; Duan et al., 2006). A review by Chua et al. (2008) 
reports that in healthy aging, DTI abnormalities occur in the frontal re- 
gions, specifically the frontal white matter, anterior cingulum and the 
genu of the corpus callosum, while in AD, DTI abnormalities are concen- 
trated in the posterior regions. Some studies have dissociated the region- 
al effects of age and dementia, with age effects greater in the anterior 
corpus callosum and frontal white matter, suggesting an anterior-to- 
posterior gradient, while individuals with early-stage dementia exhibit 
minimal additional frontal deficits, but rather show greater white matter 
deterioration in posterior lobar regions (Head et al., 2004). 

As in prior studies of aging and AD (Sullivan et al., 2010; Acosta- 
Cabronero et al., 2010), FA was the least sensitive measure when com- 
paring diagnostic groups and cognitive score associations in both ROI 
and VBA analyses (Fig. 2). A recent study found that FA measured 
using TBSS was more strongly associated with spectroscopic measure- 
ments than FA measured using voxel-wise averaging which is suscepti- 
ble to partial-averaging artifacts (Wijtenburg et al., 2012). However, we 
found that TBSS ROI results were not substantially different from 
data computed from averaging diffusion indices in full ROIs (Inline 



Supplementary Fig. S2). A recent AD DTI review - based on 55 different 
studies - noted that MD values have more discriminative power than FA 
values and higher effect sizes for case-control differences in the frontal, 
parietal, occipital and temporal lobes (Clerx et al., 2012). To explain age 
related increases in MD without significant FA changes, Zhang et al. 
(2011) suggested that brain degeneration in aging may be caused, in 
part, by tissue damage due to processes such as focal ischemia. This 
may result in lower tissue density, increasing water diffusivity but 
maintaining underlying directional structure. 

FA and MD are summary measures based on the ratio and mean of 
the eigenvalues respectively, but AxD may reflect axonal injury, and 
RD may reflect demyelination (Song et al., 2003; Song et al., 2005; Sun 
et al., 2006; Hofling et al., 2009). Here, when comparing NC to AD pa- 
tients in both ROI and VBA analyses, RD was the most strongly associat- 
ed with WM deficits (Fig. 2a-b), followed by MD. Increased RD in AD 
relative to controls, therefore, may reflect demyelination in AD. As in 
prior studies that reported AxD increases in normal aging and 
Alzheimer's disease (Fellgiebel et al., 2004; Sullivan et al., 2010; 
Acosta-Cabronero et al., 2010; Agosta et al., 201 1 ), we also found higher 
AxD in MCI, AD, and associated with clinical impairment. In fact, we 
found slightly larger effect sizes for AxD when comparing CN and 1- 
MCI subjects in the VBA analyses, and CN and MCI subjects in the full 
ROI analyses (Table 3b). Recent longitudinal and cross-sectional studies 
suggest that AxD might be more sensitive to detecting early changes, 
while RD becomes progressively better as disease progresses 
(ODwyer et al., 2011; Acosta-Cabronero et al., 2012). Ultimately, the 
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Table 4 

Atlas ROI average anisotropy and diffusivity. Here we highlight regions with the 10 lowest FDR significant (Benjamini and Hochberg, 1995) p-values (greatest effect sizes) when assessing 
cognitive test score associations in the entire study population. We also note the total number of significant ROIs (after FDR correction) out of the 43 tested and the critical FDR p-value. 
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Atlas ROI average anisotropy and diffusivity. Here we highlight regions with the "top 10" FDR significant (Benjamini and Hochberg. 1995) p-values (greatest effect sizes) when assessing 


cognitive test score associations within diagnostic subgroups. We also note the total number of significant ROIs (after FDR correction) out of the 43 tested and the critical FDR p-value. 
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Fig. 6. Average FA and RD Y axis looks so much bigger than Average AxD and RD Y axis of the 43 ROIs, these 9 ROIs were consistently more sensitive to detecting differences across dif- 
fusivity and anisotropy measures in at least two cognitive and/or diagnostic analyses. These graphs show the average FA, MD, RD and AxD in the ROIs in each diagnostic group. "TOTAL" 
refers to the ROI generated by combining all the atlas ROIs. Overall, we see a decrease in anisotropy and increase in diffusivity with each stage of the disease. 



the study design. Additionally we did not adjust for WMH burden, or ex- 
plore the extent the DTI metrics corresponded with WMH. Finally, the 
ADNI data set is constantly re-evaluating subject diagnosis, cognitive, 
and basic demographic data. The data evaluated here was what was 
available at the time of download and may be modified by ADNI. 

5. Conclusion 

DTI offers an extensive set of biomarkers for disease detection and 
monitoring of cognitive decline. We found anisotropy and even more 
widespread diffusivity disruptions in projection, association, and com- 
missural tracts, particularly in tracts that pass through the temporal 



lobe and posterior brain regions (especially the left hippocampal cingu- 
lum), in elderly AD and MCI patients. These disruptions were associated 
with neuropsychological and cognitive deficits. As ADNI2 progresses, 
new subjects are scanned and new measures of WMH are being added 
to the database. These results should be verified with larger sample 
sizes and the relationship between DTI based WM disturbances and 
WMH can be further resolved. Future tract-based connectivity studies 
may also shed light on how integrity in the tracts affects regional gray 
matter. We may also be able to create a statistical region-of-interest 
that may outperform atlas based ROIs (Hua et al., 2013; Gutman et al., 
2013). ADNI2 is a longitudinal study that will eventually allow us to in- 
vestigate which of these subjects develop AD, and if these early WM 



Table 6 

Cognitive correlates of the average anisotropy and diffusivity in different atlas regions of interest. Here we highlight regions with the "top 10" FDR significant (Benjamini and Hochberg, 
1 995) p-values (greatest effect sizes) when assessing executive function (ADN1-EF) and memory (ADNI-MEM) composite neuropsychological score associations in the entire study pop- 
ulation. We also note the total number of significant ROIs (after FDR correction) out of the 43 tested and the critical FDR p-value. 





FA 




MD 




RD 




AxD 




ROI 


p-value 


ROI 


p-value 


ROI 


p-value 


ROI 


p-value 


a) ADNI-EF 


SSL 


5.77E-5 


SSL 


1.85E-9 


SSL 


9.15E-10 


SSL 


4.21 E- 8 




TOTAL 


1.38E-4 


SSR 


8.12E-9 


SSR 


8.81 E- 9 


SSR 


4.78E-8 




CCC 


2.46E-4 


CGCL 


4.68E-8 


CGCL 


1.39E-7 


CGCL 


1.55E-6 




SSR 


3.42E-4 


CGCR 


3.1 0E- 7 


CGHL 


3.29E-7 


CGCR 


1.75E-6 




PTRR 


4.33E-4 


CGHL 


6.04E-7 


CGCR 


1.62E-6 


RLICL 


2.17E-6 




ACRL 


4.66E-4 


RLICL 


1.46E-6 


TOTAL 


1.63E-6 


TOTAL 


4.02E-6 




CC 


4.83E-4 


TOTAL 


1.61E-6 


IFOR 


7.14E-6 


CGHL 


5.08E-6 




Fx/STL 


7.55E-4 


IFOR 


3.31E-6 


ACRL 


9.49E-6 


ECL 


5.85E-6 




PTRL 


7.90E-4 


SLFL 


4.27E-6 


ACRR 


1.08E-5 


see 


7.35E-6 




see 


8.84E-4 


ECL 


8.76E-6 


SLFL 


1.85E-5 


SLFL 


1.21E-5 


# sig/43 


28 (p<2.39E-2) 




43 (p < 3.24E- 


2) 


43 (p < 3.27E- 


-2) 


41 (p < 3.49E- 


2) 



b) ADNI-MEM 



# sig/43 







CGHL 


1.38E-9 


CGHL 


1.68E-9 


CGHL 


5.28E-9 






SSL 


4.70E-8 


SSL 


1.26E-7 


SSL 


7.59E-8 






CGCL 


4.14E-7 


IFOR 


2.41 E- 6 


CGCL 


2.88E-7 






IFOR 


5.05E-7 


CGCL 


4.24E-6 


SSR 


4.83E-7 






IFOL 


1.03E-6 


IFOL 


8.30E-6 


ECL 


8.94E-7 






SSR 


2.21E-6 


SSR 


9.67E-6 


IFOR 


1.40E-6 






CGCR 


4.87E-6 


CGCR 


3.68E-5 


IFOL 


1.44E-6 






ECL 


2.14E-5 


ECL 


2.22E-4 


CGCR 


3.15E-6 






CGHR 


1.33E-4 


CGHR 


3.36E-4 


UNCL 


4.44E-5 






UNCL 


1.65E-4 


UNCL 


3.69E-4 


CGHR 


5.00E-5 


0 




35 (p < 3.06E- 


-2) 


33 (p < 335E- 


-2) 


35(p<3.51E- 


2) 



T.M. Nir et al I Neurolmage: Clinical 3 (2013) 180-195 



193 



aberrations help predict future deficits and conversion to AD. Future 
studies combining machine learning methods with different 
modalities - including CSF and proteomic markers - may ultimately de- 
termine the best way to distinguish diagnostic groups. 
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